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Abstract 

Entropy production for a system outside the thermodynamic limit is formulated using Hill's 
nanothermodynamics, in which a macroscopic ensemble of such systems is considered. The external 
influence of the environment on the average nanosystem is connected to irreversible work with an 
explicit formula based on the Jarzynski equality. The entropy production retains its usual form as 
a sum of products of fluxes and forces and Onsager's symmetry principle is proven to hold for the 
average nanosystem, if it is assumed to be valid for the macroscopic ensemble, by two methods. 
The first one provides expressions that relate the coefficients of the two systems. The second gives 
a general condition for a system under an external force to preserve Onsager's symmetry. 

PACS numbers: 05.70.-a 
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Non-equilibrium thermodynamics tries to understand macroscopic systems out of equilib- 
rium, and particularly those in steady states, without having to resort to the (unattainable) 
dynamic description of all their microscopic degrees of freedom, using instead the same few 
macroscopic variables as in the thermostatic case but allowing for situations where they 
depend on time. Interest on this field can be traced back to the works of Boltzmann on 
Thomson's hypothesis about the division of a general process into a reversible and an irre- 
versihle part Q. However, it was not nnti. 1931 that a truly systematic derivation of the 
thermodynamics of irreversible processes near equilibrium was attained by Onsager |2J, yj 
and later refined by Casimir 4j. 

Entropy production is perhaps the most important concept in non-equilibrium thermody- 
namics, totally absent from thermostatics. It is usual to look at it as a function of two sets 
of variables, the thermodynamic fluxes {<fti} and forces {Fi}, defined in such a way that this 
production can be expressed as a sum of products of conjugates, AS = ^ F^, the fluxes 

i 

being zero at equilibrium. This expression is supplemented by a set of phenomenological 
relations which gives the fluxes as functions of the forces, these relations being such that 
the forces cancel at equilibrium. It is an experimental fact that there exists a neighbour- 
hood of equilibrium where the relations between the two sets of variables are linear, that is, 
4>i — X] LijFj. 

3 n n 

Onsager's main result [2j, [3] is the symmetry of the phenomenological coefficients Ly = 
Lji, proven on the basis of two general hypothesis: regression of fluctuations and microscopic 
dynamic reversibility. Systems under the effect of external magnetic fields or Coriolis forces 
are exceptions already known to Onsager and later treated by Casimir in Ref. |J. More 
recently, it has been shown that the second hypothesis can be dropped for certain models 

In the last decades, interest in nanoscopic systems has led to put them in the front line of 
science and technology. Important research from the point of view of statistical mechanics 
las been done during the last decade, leading to such notable results as the Jarzynski equality 

fluctuation theorems, which have been experimentally 



6| and the Evans |7| and Crooks 
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Nevertheless, since the 



verified |9|]. A good overview of these topics can be found in Ref. 
seminal work of Hill ll] in the early 60 's, not much theoretical efforts were dedicated to the 
strictly thermodynamic formalism in the nanoscale up to the first years of the XXI century, 



when the same author revisited his own work and renamed it as "nanothermodynamics" 
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13] . Hill generalized the equations for open systems introducing a term associated 



to the number of small systems in a macroscopic ensemble of them that explicitly takes 
into account the energetic contribution of surface and edge effects, system rotation and 
translation, etc., usually negligible for macroscopic systems. 

In this context, it is of great importance to have a theoretical framework for describing the 
operation of nanomachines. Therefore, some thermodynamic results have been tentatively 
extented to systems far from the thermodynamic limit. Particularly, regarding the theory of 
non-equilibrium fluctuations (of interest for the development of nanomotors) , the validity of 
Onsager's reciprocal relations is sometimes taken for granted [141 ] . It is thus desirable to put 
the non-equilibrium thermodynamics of small systems on firm theoretical foundations. The 
formulation of such a non-equilibrium nanothermodynamics is the main aim in this report. 
By analogy with Hill's equilibrium theory, the number of nanosystems in an ensemble, which 
can be modified by production, destruction and transport, is introduced as a macroscopic 
variable that survives in the nanoscopic description. 

In order to study the thermostatics of a small system, Hill H] started with a large 
number TV of them, so that the ensemble itself was a system in the thermodynamic limit. 
The author developed a theory suitable for measurement devices that interact with many of 
the small systems in such a way that the relevant thermodynamic quantities are not those 
of an individual nanosystem, but their averages over a significant number of them, which in 
a homogeneous system will be equal to their average over all of the nanosystems. In other 
words, if each nanosystem is described by a set of extensive variables {X^} _ v and the 
state of the total system is characterized by jxf^ j U {A/"}, the quantities accessible to 

( ~< \ X {a) ~\ v 

measurement are IX 1 -"- 1 := ~~fcr f • The entropy of a nanosystem can likewise be defined 
as S = j^. 

The total system obeys the usual set of thermostatic relations, particularly the Gibbs 
and Euler equations in entropic form, dS t = V^dXf 1 — and S t — X] V^Xt ~ ^r> 

a a 

with : = ( ^prj ) {/3# ) an d £ '■= — ^ (fff)x (Q) ' Obviously, in the description of the 
overall system e is simply the chemical potential associated with the number of nanosystems. 
However, in the thermodynamics of small systems it is called the subdivision potential, a new 
variable with no analogue in conventional thermodynamics. 

To formulate the Euler equation for the nanosystems it is enough to divide both terms in 
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the Euler equation by J\f. Taking into account that dS = jj [dS t — SdAf] , the Gibbs equation 
can also be formulated and, subtracting the two equations for dS, an inhomogeneous pseudo- 
Gibbs-Duhem equation arises. Thus, the thermodynamic equations for the average small 
system are: 

s = Y.y (a) x (a) -f (la) 

a 

ds = J2y {a)d x {a) ( lb ) 

a 

-<^) = -E^ ( " ( lc ) 

a 

Comparing equations (1Tb]) and (jTal) . it becomes apparent that S does not satisfy Euler's 
theorem and thus it is not an homogeneous function of {XW} in the nanothermodynamic 
formalism. The thermostatics of a small system depends on its environment through e. 
Thus, a small system has more degrees of freedom than its large counterpart. The addi- 
tional contribution to the entropy (or, equivalently, to the internal energy) comes from the 
aforementioned interface, edge, rotation and traslation effects, which must become negligible 
as the size of the system is increased, if conventional thermodynamics is to be recovered. 

Recently, Ben-Amotz and Honig [15| have used the Jarzynski equality to give a general 
expression, dS = ^ dW ^ xW + k B log ^exp , for the entropy production of a system 

under a time- dependent constraint x {t) i n contact with a thermostat at the (possibly also 
time-dependent) temperature T, averaged over the processes compatible with that constrain 
(dW is the elementary work associated to a particular process). If it is assumed that this 
kind of operation amounts to an average over the ensemble of nanosystems (a reasonable 
hypothesis since Af is large) it is possible, using (JTJ), to give an expression for the change in 
e between times and to during the process determined by x (*)> suitable for measurement 
or simulation: 

A^) = ^E^ (Q) lo° 

a 

-aw 



| ( ( nr} + ^rio g {,-x,)( TV ) ) . (2) 

x(t) 

The central part of this report is devoted to the application of Hill's course of reasoning 
to a system out of equilibrium in the thermodynamic branch (linear regime). For simplicity, 
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only the case with two homogeneous macroscopic subsystems (A and B) will be considered. 
The results can be readily generalized to an arbitrary number of partitions or even to a 
continuous distribution, as long as large enough macroscopic differential volumes are taken 
in order to assure that they contain sufficient numbers of nanosystems. The method used 
is valid as long as differential calculus can describe the changes in the variables of the 
nanosystems to a good approximation (i.e. they are not too small). 

Suppose that the systems are separated by a diathermic, permeable and deformable wall 
and slightly out of equilibrium with each other. With the total system AU B completely 



Ma + Mb, however, can vary since it is perfectly conceivable that the nanosystems (e.g. 
micelles) could split or merge even in a macroscopically isolated system. Conservation of 
M would be a reasonable assumption in two opposite limits: static nanosystems whose 
dissociation energy is so high that interactions with their environment cannot split them, 
and highly dynamic nanosystems which are continuously reorganizing, but in such a way 
that the fluctuations in the total number of systems are small compared to the average 
value. It is possible to choose the time derivatives of all extensive thermodynamic variables 
as fluxes. Expanding AS t — S t — S tjeg = AS At + ASst to second order in these coordinates: 



where AX denotes the deviation of X from its equilibrium value and all the derivatives are 
evaluated at equilibrium. Note that the bars have been dropped for notational simplicity. 
This implies that y^ = y^ for all a and EaAMa — —£b^Mb and, given the conservation 
of X^ a \ it follows that AX^ = -AXgl Taking this into account and differentiating the 



isolated, the total deformation variables 




a=l 



are conserved. M = 




(3) 
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previous expression with respect to time: 

~dy [a) 




AAf A 



dy B 



(a) 



dAf, 



-AAf B + 



+ 



dy { B ] 



dX 



iff) 

Bt 



B 



AX 



(3) 
At 



AXif + ^AAf A 



(a) 



dVB AX 



^ dAf, 



B 



(a) , 



6W| 



AAOi 



AAf A + 

axt b . 



(4) 



This expression has the form of a sum of products of fluxes and forces, the terms inside 
square brackets being the forces and F^, To translate this expression into the 

nanoscopic language, the following equalities must be used: 



(a) 



AX 



(a) 
It 



AS 



dA (M I X ( I 
dt 



AS t _ AAf A 
>->t _ 



AX 



(a) 
it 



A/" 1 XT 2 
AfjAX^ + zj a) AM - AXi Q) AA/}, 



(5a) 

(5b) 
(5c) 



giving the result 



AS = F (a) AX ( f + F Ma AN a + F Ub AU b , 



(6a) 



with 



* Af 



F^t + EF^X^ 



Fn b 



Fat, 



Bt 



AT 



Af 2 



(6b) 

(6c) 
(6d) 



Af Af 2 '' 

which means that the entropy production of the average small system can also be written as a 
sum of products of fluxes and forces. As mentioned previously, there exists a neighbourhood 
of equilibrium in which a set of linear phenomenological relations between these variables 
holds. By means of equations (jHJ) and (jEJ), the macro and nanoscopic linear coefficients can 
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be related: 

rjw = ^ [m a l^ + xfL^) + 

+ ^ aL («Vx) + jtf) £ (V^) ( 7a) 
= J_ (tf A LW + X A a) L^'^ (7b) 
4^°) = _L (.A/^ (A/>a) + jf^M^)) (7c) 

(7d) 

The macroscopic system satisfies Onsager's reciprocity by hypothesis, i.e. L^f^^ = L^ Aa \ 
It is easy to see, starting with the last equations of the previous block and progressively back- 
substituting, that the nanoscopic coefficients are also symmetric in this situation. These 
proportionality relations ensure that the second law of thermodynamics is obeyed by the 
average systems (although it can be transitorily violated by a small system) a topic also 



discused in Ref. 
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The quantity St is indeterminate in one additive constant; therefore, the component of 
the fluxes proportional to it, arising from the terms in ( l6c|) and ( I6dl) . must be zero. This 
last condition is equivalent to: 

These equalities allow for further interpretation of (J7J) using the change of variables 
(Ma, Mb) — > (A/", D), with D := Ma — Mb- The time derivative of M represents the creation 
of nanosystems per unit time, while at fixed M the time derivative of D corresponds to 
the transport of small systems. The two possible causes of variation of (Ma, Mb) are thus 
decoupled by this change. Furthermore, application of the Curie principle shows that vari- 
ation of D can only be coupled to vectorial fluxes such as those treated in this report, while 
variation of M can be coupled with chemical (scalar) processes. In particular, they cannot 
be coupled with each other, because of the different tensor rank of the forces involved. 

Equations (JSJ) are equivalent to stating that L (aA ° = L {olMa \ L^ dd ^ = L^ MaMa ^ and the 
rest of the phenomenological coefficients involving M or D are zero. L[ DD ^ is proportional, 
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with a factor to its nanometric equivalent. The remaining coefficients are: 

T (M A a) , T {M B a) 0K r 
L (ff«) = Lt +Lt = 2A(_ LCVa) ( } 

T (MAT) Y (a) 

LfD) = (9b) 

The last equation can be used as a definition of in the particular case in which 

Af is fixed. This result means that the contribution of a flux of nanosystems to that of 
an extensive variable X^ is proportional to the amount Xj£ of that variable that each 
nanosystem carries with it in its transit from A to B, as expected. 

The above results can be reproduced by an alternate method closer to the well- 
known macroscopic proof. This method rests on the hypothesis of regression of fluctu- 
ations, which states that the reaction of the system to a small deviation from equilib- 
rium caused by an external force is the same as if it was caused by a spontaneous fluc- 
tuation. Moreover, if {-Z^}"^, := jx^i U {Af}, dynamic reversibility can be ex- 

a ~ L J a=l 

pressed as (t) AZ™ (t + r)) = (AZ^ (t + r) AZ^ (t)) , which straightforwardly 

gives (^AZ^ AZ^^ = ^ AZ^AZ^j. Substituting here a linear development for a system 

in the thermodynamic limit, analogous to that of previous sections, AZ^ a = ^ijf^z^, 

P 

with := ^py, and using the fact t 
macroscopic system, see for instance Ref 



.a 



rat ^( a 'z^) = —IcbS^' (easy to prove for a 
\) the symmetry of the phenomenological matrix 
follows immediately. However, for the nanosystem, the temporal evolution of the internal 
variables cannot be related only to their deviations from equilibrium, adopting instead the 



more general form 161 ] : 



Al(Q) = E iM aS)+/ (a) W- ( 10 ) 

p 

where (£) represents a general force acting on the system. From Eq. (I5c|) . AX^ = 
— - — j^ t,eq — . Approximating ~ and taking into account that derivatives at constant 
Af are equal for the macroscopic and average system, it follows directly that (Al^^py) = 
~ kB M ^ ■ Thus, temporal reversibility implies that: 

L {aP) _ L (M = (f(a) AX (P) _ f(P) AX ^) . (11) 

This result represents a general condition for Onsager's symmetry to hold when evolution 
of the system is conditioned by a general external force. In this particular CctSC, clS it has 
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previously been shown, the external forces are proportional to AAf through constants 
£(oA0 ; which makes the right hand side of the previous equation trivially equal to zero, 
thus proving the symmetry of the phenomenological submatrix involving only the internal 
coordinates. This second method cannot prove anything about the coefficients involving 
Af, an external variable. To show their symmetry it is still necessary to relate them to 
their macroscopical counterparts. This formulation, however, connects more easily with the 
language of the fluctuation theorems. 

New cross-transport phenomena associated with the new degrees of freedom must appear 
in the nanosystems. As predicted by Eqs. (Q, if the number of small systems is kept con- 
stant, these phenomena will consist simply in the exchange of extensive variables transported 
along with the nanosystems. Dynamic nanosystems, such as micelles, could be thus wor- 
thier of study. Since Hill's equilibrium formalism has already been applied (and success fully 
compared to experimental data) to nanostructures such as nanosolids and nanowires 17| . 
they could also be good candidates to find these cross-phenomena; for instance, in electric 
or thermal measurements. 
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